function H = gaussianFilter(fftStruct, sigma)
% H: Filter transfer function
% sigma: std deviation

[H r c] = initFilterSize(fftStruct.fftImage);

%squared std dev
sigma2 = sigma*sigma;

for i = 1 : size(H,1)
    for j = 1 : size(H,2)
        dist = sqrt((i - (r/2))^2 + (j - (c/2))^2);
        H(i,j) = exp(-(dist^2)/(2*sigma2));
    end
end
end